NASA Contractor Report 201584 
1CASE Report No. 96-45 



MAXIMUM PRINCIPLES AND APPLICATION 
TO THE ANALYSTS OF AN EXPLICIT TIME 
MARCHING ALGORITHM 


Patrick Le Tallec 
Monday D. Tidriri 


L 


STATEMENT & 

Ass novae tot puroic ifi&ojud 
Otummctt Usuuxmtd . .... 


NASA Contract No. NAS1-19480 
July 1996 

Institute for Computer Applications in. Science and Engineering 
NASA Langley Research Center 
Hampton, VA 23681-0001 

Operated by Universities Space Research. Association 



National Aewnautics and 
Space Administration 

Langley Research Center 
Hampton, Virginia 23681-0001 


DTIC quality inspects© a 

19960924 141 


THIS DOCUMENT IS BEST 
QUALITY AVAILABLE. THE COPY 
FURNISHED TO DTIC CONTAINED 
A SIGNIFICANT NUMBER OF 
PAGES WHICH DO NOT 
REPRODUCE LEGIBLY 



Maximum Principles and Application to the 
Analysis of an Explicit Time Marching Algorithm 

Patrick Le Tallec 
INRIA 

Doraaine de Voluceau Rocquencourt 
B.P. 105 Le Chesnay Cedex France 

Moulay D. Tidriri * 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton VA 23681-0001 


Abstract 


In this paper we develop local and global estimates for the solution of 
convection-diffusion problems. We then study the convergence properties of 
a Time Marching Algorithm solving Advection-Diffusion problems on two 
domains using incompatible discretizations. This study is based on a De- 
Giorgi-Nash maximum principle. 
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1 Introduction 


Domain decomposition methods have recently become an efficient strategy 
for solving large scale problems on parallel computers ([1], [2], [3], [4], [5], 
[6]). Nevertheless, they can also be used in order to couple different models 
jll], [18]. In this paper we will examine a domain decomposition strategy 
which can be applied to such situations. 

This approach was introduced in order to solve several difficulties that 
occur in fluid mechanics. In particular, our aim is to introduce several sub- 
domains in order to do one of the following ; 

• Solve different problems on each subdomain. 

• Use different kinds of approximation methods on each subdomain [7]. 

• Use “local refinement techniques” or “mesh adaptive techniques”, lo- 
cally, per subdomain ([10]). 

The subdomains fully overlap and the coupling is achieved through “fric- 
tion” forces acting on the internal boundary of the domain, these friction 
forces being updated by an explicit time marching algorithm. 

Several versions of this methodology have been studied in [15]. In [15] the 
emphasiz was on the implicit time discretization version of this algorithm, 
we focus in this paper on the explicit version of this methodology. The 
theoretical study of our method will be done on an Ad vection- Diffusion 
problem, which will serve as our model problem. The analysis will be made 
at the continuous level, independently of any (space) discretization strategy, 
which means that the derived results will be mesh independent. 

In the next section we develop a maximum principle for general second- 
order elliptic problem based on the De-Giorgi-Nash theory. In section 2 and 
3, we develop estimates for the solution of the convection-diffusion problems 
respectively with DiredJet- Neumann and Direchlet boundary conditions. 
These results are based on the maximum principle of section 2. We then 
apply these tools to the analysh of an explicit time marching algorithm. We 
also study a fixed point method for the implicit time marching algorithm of 
[15]. Practical applications of the lime marching algorithm to real life CFD 
problems can be found in [14], [19], [20], and [21]. 
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2 Local estimates 

In this section we shall establish a maximum principle for an arbitrary elliptic 
operator of second order. These tools are central to the development of 
our theory in order to derive the convergence analysis of the explicit time 
marching algorithm described in section (5.1). 

Let L be an operator written tinder the form 

Lu = a**(z)DijU + b'(x)DiU + c(x)u , 

for any u in W 2 ’ n (ft) ; with ft a bounded domain of IT. The coefficients 
a?Kb x and c,i,j = 1, n are defined on ft. As usual, the repeated indices 
indicate a summation from 1 to n. 

We suppose that the operator L is strictly elliptic in ft in the sense that 
the matrix A of coefficients is strictly positive everywhere in ft. Let A 
and A denote respectively the smallest and the largest eigenvalue of A. Let 
V denote the determinant of the matrix A and V* = X) 1 /". We have 

0 < A < V m < A. 

We suppose in addition that the coefficients a ij ,b' and c are bounded in ft. 
and that there exists two positive real numbers 7 and £ such that : 


A/A < 7, (L is uniformly elliptic) (1) 

(l*l/A) a < (2 ) 

Now, we are in a position to state; the principal result of this section, 
proved in annex. 

Theorem 2.1 Let u <S W 2,n (fi) and suppose that Lu > / with f € L n (Cl) 
and c < 0. Then for all spheres B -■ Bw(y) of center y and radius 2 R 
included in ft and for all p > 0, we have : 

S CW(]5j / B (« + )') 1/p + f II/II„.b}, (3) 

where the constant Cr depends on (n,y,SR 2 ,p), but is independent of c. 
Above u + = max(u,Q). 
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Remark 2.1 The statement of the same theorem can be found in [12], un- 
der the assumption 


jc|/A < 6. (4) 

So, there the constant Cr depends indirectly on c through 6 - That is 
exactly what we want to avoid , since we would like this constant to be inde- 
pendent of c (see section 5. 1 ). 


3 First fundamental estimate 


Let ilioc be a connected domain of Bl n with Sli oc C O (Figure 1). The 
boundaries of the two subdomains are defined as follows: 


T& = dft n dClioc, ( interntd boundary) 

r, = dilioc H SI, ( interface) 

Toe, - (farfield boundary) 

We denote by n the external unit normal vector to dSl or dSli 0C . 

• et V be a given velocity field of an inviscid incompressible flow such that: 


divV as 0 in SI, 
V.n — 0 on Tj. 


(5) 


We shall derive an estimate for the solution of the following Direchlet- 
Neumann problem: 


Cv = -i/Av -f V ■ Vv + in SI, 

T 

V = 0 on Too, 

dv 

SZ ~ S on T b , 


( 6 ) 

(7) 

( 8 ) 


where the function g is given in J?" 1 ^^) and the coefficient r is strictly 
positive, and u is the diffusion coefficient. Let W be the sub-spa.ee of ff 1 (Sl ) 
defined by 


W = {tu 6 H'{Sl) | u> = 0 on F<*,} 


( 9 ) 
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We then define the following bilinear forms on W 



v'VvVw + 



( 10 ) 


(«,«>) = J vw. ( 11 ) 

The first basic problem associated to (8)-(8), can be written as follows: Find 
v € W satisfying : 


a(v,-w) + (l/r)(v,w) — j gwdT, Vtn 6 W, (12) 

Moreover, we assume that the coefficients v and r satisfy the following 
relation; 

vr<l. ( 13 ) 

This hypothesis is not necessary but simplifies the proofs to come. More- 
over, it is not restrictive, since we would like the convergence for small t 
( see section 5.1). 

L«*t d denotes the overlapping distance as described in the Figure 2. Let 
then 0 be a real number such that 


0 < 0 < 3 \fvjd, 

and set 

k = ft /(uy/P). 

The first basic result states the global H l estimate of the solution of the 
first basic problem (12) in terms of the boundary data g. 

Lemma 3.1 There exists a constant c„ such that we have: 

IM!l,o < (e«>A')||<7||-.i/2,r k , (14) 

Proof of lemma. 3.1 


By using the relation (5) we have the following equality: 


udiv(Vv) = 1/2 J div(VV) 



= 0, V v € W. 
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(15) 


Choosing w = t; in (12), we til'**! obtain 

/ {i'jVu| 2 + (l/r)v 2 ) = / gv. 

J n 

From this equality we deduce the following estimate: 

^ IL?l)-i/2.rJl w fli/2,iv 
The application of the trace theorem yields the estimate (14), which implies 
in particular 

IMko < {colv)\\s\\-H 2 X h . (16) 

■ 

Let Slj be the subdomain of width ~ with external boundary T,- as described 
in the Figuie 2. Let K y ~ Bd(y) be the sphere of center y and radius |. 
There exist yi, . . ■> yj belonging to fi, such that 

: = C Uj = i A' yj . 

We define then K by setting 

K = U I'.,*-,,,. 

Tbe next lemma states the local estimate of the solution v of the first basic 
problem (12). 

Lemma 3.2 There exists a constant ci such that: 

IMU.fif < Ci||l>||o,(i. (17) 

where Ci is a constant depending only on }',i,6d 2 and (3/2 d) n ^. 

Proof of lemma 3.2 
The operator 

L = -C 

satisfies the assumptions of the theorem 2.1, with c = -1 jr and / = 0. 
Applying then this theorem with p = 2, y € £1, we obtain 

IMUo.ffy < CllM!o,£ i(//4 (y)- 


5 



Therefore 


( 18 ) 


IMIoo.ft'j, < Cj IMjo.fi, 

where cj is a constant depending only on v, 7, id 2 and (3/2 d) n /2. 

Applying the relation (18) to each K V) we obtain: 

IMU.K < sup Ci J ||v||o,fj. 

Setting ci -- sup c ljs we finally have 

i=i 

IMU,*" < c ? ||nj|o,n- (19) 

And the lemma is proved. 

■ 

We shall now establish other local estimates for the solution v of the first 
second basic problem. For any Mi in fl,, we introduce (see Figure 2): 

« Bi = the ball centered on M, of radius d/ 6 , 

• vi = <®p[b(r 2 - d 2 /36)j||v|J 00i aj Bi . 

We then have: 

Lemma 3.3 The solution v of the first basic problem satisfies: 

l v (Mi ) | < exp(-fcd 2 /36)||o|| 00 ,a B ., VM,- € H t . (20) 

Proof of lemma 3.3 

The operator C applied to v,* can be written in polar coordinates (with 
r ^ M X M) : 

lv, ss 4(-fc 2 */r 2 - kv + ^v.e r r 4- -i-W. 

z 4r 

Therefore 

iv, > 4(-* 2 i/r 2 - “|V.e r |r 4 (~ - fcr/))v,. (21) 

We set then : 

9(r, k) ~ a(k)r 2 -f 6(b)r 4 c(&), (22) 

with 
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a(k) = —k 2 v 
m = -||V.e r j 

c(fc) - -p — kv. 

At 

We seek to satisfy the following relation : 

d 

C < inf <p(r , k ) for 0 < r < 

o 

As <p(r., k) decreases on 1R + , this will be satisfied iff 

<p(d/6) > 0, 

i.e. iff 


fcW kd\\V\\ .1 

36 1 — 

We replace k by its value. Therefore, we have to satisfy 

W _ /?rf||V|| , J_ „ J*V_ > ft 
(36z/r) 12.ty7 4r i/a/t “ 

Multiplying by </r, it follows that 

4V? 1 9«- ’ S M 12t* '• 

The constraint /? < Z\fvjd, finally yields after division 

V (r, *) > 0 iff > />[! + ][1 - ?£)-' ■ (23) 

From the relation (21) and the previous calculation, we deduce that for 
/? < 3 \fvjA and r satisfying the above inequality, we have 


Lvi > 0 = Lv. 
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In addition, by construction 


t>{ > v on dB{. 

Consequently, by using the maximum principle we obtain the following re 
lation : 


v < Vi in 2?,. 


In particular 

v(Mi) < cx-p(-kd 7 

We do the same for —v> and finally we have 


K-^i)! < exp(-£d 2 /36)|M| 0i ,as,, VAf, e ft,-. 


(24) 


Let fioo be defined by Hoc- — H \ fhoc The next result establishes an H * 
estimate of the solution v of the first basic problem on the domain 

Lemma 3,4 There exists a constant r.$ #uch that: 


Proof of lemma 3.4 
Let £ € J9 rl (fl) be such that : 

/ £=1 in fU, 

1 supp^ c fi. Ufloo. 

We have using (12): 

I (-vAv + div(Vt>) + v/t)£ 7 v = 0. 

J Q 

By using the Green's formula we deduce : 

/ -vAvfv == / HV(£»)| 2 ~ / r|V<f| V. 

Jo JQ 


(25) 


(26) 


( 27 ) 
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On the other hand, we have : 

/ div(V?,-)f 2 v = / dhWV/2) - / y.V^ 2 . (28) 

J Cl Cl J Cl 

Using the relations (27) and (28), (26) becomes 

0 = j[HV(^)| 2 + div(U^ V/2) -f fvt/r)- j (w 2 \Vtf + v 7 tV-VZ) 
= Kl v (^ v )l 2 + IH 2 ) + J n ^ l ! r - v^v 2 - jf (*a* a |V£ 2 + v 2 cf.vo 
= X *<^ + W 4 ) + X.K|V(^)| 2 -HH 2 )4 foVT-uW 
- I (v» 2 \Vt\ 2 + v*tV.VO- 

J fit 

Hence, we obtain : 

"IMIi.noo + X z/ (! v ^ v )! 2 + l£«l 2 ) + X^/ T _ ^^ 2t;2 = 

f (^IVCI 2 -f v 2 {V-V4). 

J Oj 

The relation (13) then yields 

< ( (™ 2 |V{| 2 + v z (V.V() (29) 

JQ % 

< !M&.q, I MV<£| 2 -KKVO (30) 

< IMILaKIW* + ||f||o < o,||V|/|{|i,n < ). (31) 

If we take £ such that 

llflkn, < 1, 
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l£li,n, = c 2 /d , 

where c 2 is a constant, (31) then becomes 

< lblU,ni\/c^ ^1 + , (32) 

which is the conclusion of our lemma. 

■ 

Now we are in a position to state the main result of this section. 

Theorem 3*1 Jjet v he the solution of the first basic y^oblevn ^ J 2J . Jf t is 
sufficiently small, we have 


(Mli/ 2 ,r, < <?i\/5/d(i +~\\V\\ aaS fdjC^ ^ 
(l/i/)erp(-fcd 2 /36)||pi|_ 1/2 , rj( , 

where C\ and C\ are constants, with C\ depending only on d,v, 7 and 6, but 
not on Ts 

Proof of theorem 3.1 

The proof of this theorem is based on the above lemmas. Since dBi C K 
We have 


iMloo.&Bi < IMU./C, (33) 

The lemma. 3.2 then implies 

Nloo.asi < Ci||u|!o.o- (34) 

Lsing the lemma 3.3 and the above estimate we obtain: 

WMi)\ < exp{—kd 7 fZ($)c\ (|v(|o,n, VJlf, € Qj. 

Consequently we have 

MU,o, < eip(-fcd 2 /36)ci IMIo.n- (35) 
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Applying the lemma, 3.1 we obtain: 


MUa < ~ezp(-*d 2 /36)|M|_. l/2JV (36) 

The application of the lemma 3.4 then yields: 

iMkfu < c^yfri (l + M^/d/c 2 ) ' (37) 

(l/i/)eip(-*d 2 / 36 )j| 5 |i_ 1 / 2 > r k . 

To conclude we use the trace theorem which yields 

IMii/2,r, < c 3 ||u||i,o eo . 

Consequently, we have the final estimate: 

!MIi/2,t\ £ c oCiCz\fci/d(\ + -~v^/c 2 ) 1/2 
(\/v)exp{-kd*[Z6)\\g\U, xri , 

which corresponds to our theorem witii C\ = c 0 C!C 3 and Cj - c 2 . 

■ 

4 Second fundamental estimate 

In this section wo shall derive an estimate of the solution of the following 
Ditechlet problem: 


-vAv + V ■ Vv + -v in fi/ oc , 

T 

(38) 

h , on l"i, 

(39) 

0, on T b , 

(40) 


where the function h is given in /3 rl / 2 (r,), the coefficient r is strictly positive, 
and v is the diffusion coefficient. The velocity field V is given by the relation 
(5). Let W be the sub-space of /7 5 (^(oe) defined by 
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W - {w € H\Skoc)\ w=:Oon r 4 }. 

We then define the following bilinear forms in W: 

o(t;, w) = v / Vv.VuH- / div(Vt>)to, (41) 

•'toloc •'ftloc 

(v,w)~ / vw, (42) 

•^O(oe 

with v and w in W. The second basic problem associated to (38)-(40) 
corresponds to the following problem: 

Find v € W such that 

a(v,vj) + (1 /t)(v,w) = i/~w, Vtt; e W , (43) 

i)|r, = h, (44) 

where: h is given in T,). We first have the following lemma ; 

Lemma 4.1 For r sufficiently small, we have 

a(w,w) + (]/r)(w,v;) > (t//2)||io||J (0toeI Vtt- € W. 


Proof of lemma 4.1: 

Under the hypothesis 1/r > vjl + (l/2/>')j|l/||^ OI and using the Ca.uchy- 
Schwarz inequality, we obtain : 

a(v,v) + (l/r)(u,t;) = j vVv.Vv + / V.Vw?) + (1 /t) / v 2 

* ^(oc ^loc v 

> v\iVv\\la + (1/t)||u||o,j - ||V|| 00 ||Vt»||o l3 ||i»||o,a 

> H|V< 2 + (l/r)||v||g, 2 -( t //2)|iVu||^ 

-(Wll^lLlMlk 

> (^/2)(||Vu||g i2 + ||v||g i3 ). 
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We will also make the simplifying assumption (.13). We first establish a 
global estimate for the solution of the second basic problem. 

Lemma 4.2 The solution v of the second basic problem (48)- (44) satisfies: 

Ml'-"- < 2( 1 + r- 2 )‘/^l + l±|pi ) ' ||fc|| 1/2 n (45) 

Proof of lemma 4.2: 

Choosing u = t; in (43) we obtain : 

V L ! Vv ! 2+ i (div(V«)® + (l /r)v*)= f v %Lh. (46) 
•'Mice Jr, on 

The lemma. 4.1 then yields 

MliAoe ^ 2 \\ dv l dn \\~\/2X,\Mn2,r,. (47) 

We shall, now establish an estimate of j|dt>/drc||_ 3 , r . Combining (43) and 
(5) we obtain: 

f dv t 1 

JT ~dn W ~ Jn ■*" (4)v')V,Vvw 4- — vw). 

Therefore, for any w in W, we have 

1 J^H ^ ll Vv llo.n lo J|V?n|!o,o Jec + (l/^)il^||oc||Vw|l 0tnioe |jtn|| 0 ,fi loe 
IMIo,n jae !MIoa« 

< (IF’lftnu. + (l/^)|iniLl|V V j|g >ftjM + OA' 2 )||<n I J 1/2 

OI Vw ll<W + HI to,rn, - (l/r 2 )|Mft 0 | J 1 '* 

< (i + + r" s ) 1 /a ||uj||j >niM . 
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The trace theorem then yields 


\\dv/ 5n||_i/ 2 , r . < (1 + r - 2 )l /2 



i + llv r H2o \ 

^ ) 


1/2 


(48) 


Combining now the relations (47) and (48) we have 

Mist** < 2(1 + t -*?* 2 (1 + ' ii^iii/ 2 ,r, m 

and hence in particular 

IWkrw < 2(1 + r" 2 ) 1 / 2 (l + 1 j| fc j, 1/wv ( 50) 

M 

Let K y — Bi/ A (y) be the sphere centered on y and of radius d/4, with y 
belonging to Tv (see Figure 3). By construction, Tv is the center surface of 
Qi oe and ft, is the subdomain of width $ centered on Tv . 

We have the following lemma: 

Lemma 4.3 There exists a constant ci svch that: 

NU.fr < ci |]v|| 0 ,n Ioc . (51) 


Proof of lemma 4.3: 

Following the same argument as in the proof of the lemma 3.2 we obtain: 

IM|oo,k„ < «i||t>)|o,n Joe , (52) 

where c\ is a constant depending only on d, and />. On the other hand 
there exist j/i , . . .,yi in ft; such that 


n* = U bm cU^ ; = k. 

V€ 0 , j * al 

By applying the relation (52) to each K Vjt we obtain 
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Figure 2: Description of the local domain f). Joc and of the splitting used in 
the maj oration of the global solution. 


IMU* ^ _f»P c iilMk«W = CiNIoAoc- (53) 

■ 

Next we shall establish another local estimate for the solution of the second 
basic problem (43)-(44). For any Mi £ fti, we introduce (see Figure 3): 

• a ball Si centered on M{ and of radius d/6, 

• the function t = exp[k(r 2 - d 2 /36)]|M|«>,aB ; - 
We then have: 

Lemma 4.4 The solution v of the second basic problem US)-(U) satisfies: 

KM) I < «*p(“^ 2 /36]||v|| oc ,, aBi . (54) 

Proof of lemma. (4.4): 


By construction of k (see the previous section), <^(r, k) is positive for all 
r € (0, d/6]. Then by following the same argument as in the proof of the 
lemma 3.3 we obtain the inequality (54). 
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m 

Let fit be the subdomain of fl/ oc described in the Figure 3. The if 1 global 
estimate of the solution of the second basic problem, is obtained in the next 
lemma. 



Theorem 4.1 For r sufficiently small, the solution v of the problem (4S)~ 
(44) satisfies: 
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||a«/9n|l-,/2A < CiJcJdl i + l + jpi 

v V* 




(1 + l/r 2 )exp(-fcd 2 /36)||A|j 1/2>r> ., (58) 

where Cj and C 2 are constants with C x depending only on d,v,i/ and S. 
Proof of theorem 4.1: 

Since dB, C K by construction, the lemmas 4.3 and 4.4 Imply: 

WUft < €xK-^ 2 /36)ci||v||o.fi <oe . (59) 

Furthermore by using the lemma 4.2 it follows: 

c,(l + l/r 2 ) , / 2 «p(-Jbi 2 /36)||A|l- ,, >r .. 

By using the lemma. 4.5 we then obtain: 

IMIi.ruuo, < 2 ^3 + j 1 

, nrrV 12 (6i) 


(1 + l/r 3 ) 1/ ' 2 ea?p(-Ad 2 /36)i|h|| 1Ar .. 
Before concluding we shall establish an estimate of the term 

||d»/0n||_j/ 2i rv 

Choosing w such that: 

w € H^Clioc), with w -- 0 on dClb n #ft,, 
and using (43) we obtain: 
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5 Convergence analysis of the explicit time march- 
ing algorithm 

Consider the following oUiptic problem: 

— 4 V.Vi i> — u&<j> — 0 in ft, 
r 

<t> = on Too, (63) 

4> =• 0 on r 6 , 

that we would like to solve by the fundamental algorithm of [15]. This 
algorithm can be written in this case as 

• set 4>° loc ~ 4> 0 i and 4>° - <j> 0 . 

• then, for n > M]‘ oc and <£ n being known, 
solve 

ai n+1 

+ KV^t 1 - vM? o y = o in ft (oe , 

^ on r„ ( 64 ) 

** 0 on F fc , 
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- t/A<f> n+1 

= 0 in ft. 


f> »+» 

= <j>oo on r^, 

(65) 

di> n+l 

- on r ‘- 


V dn 



We shall show in this section that this algorithm converges, and the 
converged solution corresponds to the solution of the initial problem (63). 
More precisely we have the following theorem. 

Theorem 5.1 Fort sufficiently small, <f> being the solution of the stationary 
problem (6$), we have : 

* \ d<l> l0 V d<f> . . 

i) —fa— converges to — - in E~ 1/2 (Ts), 

ii) </> n+1 converges to 4> in H l l*(Ti), 

iii) <^ n+1 converges to <p in 

iv) converges to <f> in II l (fl loc ). 

Proof of theorem 5.1: 

By the transformation </> n+1 with <f> the solution of the stationary 

problem, tliis problem can be reduced to the case <f>^ = 0. Multiplying the 
equation in (64) by w e W, integrating by parts, we obtain: 

L ^ * L w «‘- v «‘ 

f <W +1 

= v J Vu> e w - (66) 

We now apply the theorem 4.1 and we obtain 

II H-i/2 A < + + (67) 
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(l 4 


(1 + \/^ 2 )exp{-kd 2 jZ<o)\\4> n \\ij-j.j- t . 

Ob the other hand, multiplying the equation in (65) by t v € W and inte- 
grating by parts we obtain the equality 


/ f y,V<^ n+1 u- + v [ V<£ n+1 Vu> == v [ -]£-w (68) 

Jfi r Jo Jo Jr.s on 


with w € .ff^lt) and id = 0 on Too. Applying the theorem 3 1 to this 
problem yields: 


ll ^ n+1 lli/2,r, < ciy/c>/d(l+±\\V\\„y/dfay /2 

exp(-kd 2 /36)\\d4>l' 0 + J /dn\\_ 1/: , Xb . 

Combinig (67) and (69), we then ha.ve: 

m?oV/M\-U2T b < + j- (i + \\v\\D) 

/ 1 j \*/ 2 

fi + ~m„jd/c'A 


cisfi7/d (l + i||V| u/d/c^ 1 ' 2 

(1 4- 1 /r 2 )exp (~k^j \mjdn\U hXi , 


with k - Therefore for r sufficiently small, the coefficient of reduction 

will be dominated by the exponential term and will then be strictly less 
than 1, implying the linear convergence to zero of 
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This corresponds exactly to the statement (i). This statement combined 
with (69) leads to the convergence of <£ n+1 to 0 in ^ 1/2 (r,). Applying (14) 
with g — dtf*}/8n, we have in addition 

and therefore ||^ n+I ||i, 0 converges torero at the speed of JS^JJV&tll-t/Mv 
Applying now (45) with h = <f> n , we also have 

IK'lltA,, < 2(1 + Ur 2 )' 11 (l + ^(1 4 l|V!|^)) 1/2 ||^|| I/2 , ri . 

And then (|<^ n+1 |)x,fj also converges to zero at the speed of ||< 4 > n ||i /2 r • 

■ 


5.1 Convergence of a fixed point method for the implicit 
time marching algorithm 

The implicit time marching algorithm of [15) couples the global and the local 
problem. To uncouple them, it is advisable to use the fixed point algorithm 
below : 

• set 4>°oc,o = i'oi and <f>° ~ 1 p 0 , 

• then for k > O,^^ 1 being known, 
solve 


At 


+ div(«^+\ +1 ) - vA<t?£ k+l = 


4> n+l 


1 


n+1 




*+ 1 ~ 


0 in n toe , 

4> n k +1 on r„ 

0 on Tfc, 


(70) 


^.'+1 


-r 


At 


+ div(««0 - I'A^+I 


- 


0 in fl, 

4><x> on r w , 


-WK }/«» = on r,. 

( 21 ) 
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We will study now the algorithm (70)-(71). By setting 


4w., = - «£J+i. I7r > 

*>*.,= (« +, -< +1 ). ( 73 ) 

we have that il>i oc ,k,q and verify the following equations : 


— Z''A^j0 C| fc,g 

1! 

o 

P 

^loct 

^hc,k t q 

jj 

?r 

-i on 


0 on 

r i b, 

uAll?k,q = 

0 in H, 


tk.ff = 

0 on Too, 


,.frk.q 

'd'l’loc.k.q 

on r i. 

On 

dn 


(75) 


If A t is sufficiently small, we can apply the analysis of the previous 
section to this algorithm and we conclude that ipk t<J and tpiac conv rge 
linearly to zero. Hence the sequences 4>% +l and are Cauchy sequences, 
which converge linearly to the unique solutions <£ n+1 and of the implicit 
scheme. This guarantees tlie convergence of the above fixed point algorithm. 
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Appendix 


The main result of this section relies on the notion of a. contact set. If u 
is a continuous arbitrary function on ft, the upper contact set, denoted T + 
or T+, is the sub-set of ft, defined by 


r + = {y € ft, 3?(y) € fft n such that a(ar) < u(y ) + p ■ (x - y) Vz e ft }. 

(76) 

We see that v. is a concave function on ft iff r + - ft. When u 6 C 1 (fi) 
we must have p =■ Dn(y ) in the relation (76). In addition, when u € C 2 (ft), 
the Hessian matrix D^u = [Dijti] is negative on T + . In general, r + is closed 
in ft. 

If « is a continuous arbitrary function on ft, we define the ’'normal 
mapping” x(y) ~ Xu(y) at point y £ ft by 

Xiy) = {p € IR n , u(z) < u(y) + p.(x - y) Vz € ft}. (77) 

We can see that x(y) is non empty iff y £ T+. In addition when u € C^ft), 
we have x(y) == Du{y) on F + ; in other words y is lb' radient field of u on 

r+. 

As a particular case of the Bakelman- Alexandrov ([8] and [9]) maximum 
principle, we have under the above notation. 

Lemma .1 For u € C 2 (ft) n C°(ft), we have : 

supu < sups o« + ~4 7 ^|!a’>'D i> ti/2>-|| n(r+ 
n riWn 

with d the diameter of ft and w n the volume of a unit sphere in R n . 

For further details see [12]. 

We now proceed to the proof of Theorem 6.1, by following the steps of 
[12]. We take B = (0) and the general ease will be deduced by considering 

the coordinate transform, se — <• x - (x _ y)/2R. 

w « will begin, in first step, by showing this result for u <E o' 2 (ft) n 
W 2,n (ft) and then in a second step we will deduce the result for u € W 2 ’ n fft). 
Step 1: 
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We suppose that u 6 C 2 (Sl) 0 W 2,n (fi). For 0 > 1, we consider the cv* 
off function tj defined by 

m = (i - 


By differ enti at ion, we obtain 

D x r) ~ -2 0ii(. - |i| 2 ) 0_ \ 

Dijv = -2/My(l - !i| 2 )^ -1 + 40(0 ~ l)*.-*i(1 - \x\ 2 f~\ 
By setting 

v = nu. 

we then obtain 

a ,J bijV = rja'i DijU + 2a'i DiTjDjU + ua'i Dip] 

> T)(j - b'DiU — cu) 4- 2a^ bir/bjU 4- ua^Dijr). 


Let T+ ~ T+ be the upper contact set v, in the sphere B ; we have : 

u > 0 on T + . 

If x € 8B such that p.(x - y) < 0 we indeed have v(x) — 0. Consequently 

v(y ) 4- p-(x -- y) > v(x) - 0. 

Moreover, using the concavity of v on F + , we can estimate the following 
quantity : 

\Du\ -- (\/rj)\Dv - 

Indeed, 

|£*| < (l/»?)(|^t)|4-u|/>7j|) 

< M—pi + uN) 

< 2(1 4 - 0)rj~ i ^u. 

In tha‘ way, we have on I' - ' the following inequality : 
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-d'-’J D,jv < {(16/? 2 + 2TjQ)krj~ 7 ^Jr 

20\b\-q~^^ -f c}v 4- r}\f\. 

Since c < 0, we deduce the inequality 

-a t3 DijV < {(16,8 2 + 2 ti/3)At)~ 2/0 + (78) 

< ci\t)~ 2/0 v+ i/i, 

with ci - c(n,0, 7, b) independent of c. 

Consequently, by applying Lemma 6.1 on B, we obtain, for 0 > 2 : 

By using the relation (2), it conies 

S r * { n^ 

< c 1 d(||t ? - 2 /^|| ni6 + (l/A)||/|| n g) 

< ^411^/^11^ -h(l/A)||/|| n S ) 

< c 1 4(supu + ) 1 - 2 ^||(t 1 +) 2 ^|| n ^ + (1 A)||/|| n ,*), 

where cj is a, constant depending only on n,fi,y and 6. Here, d is the 
diameter of B(d ~ 2). 

By using the Young inequality under the form 

ab < ea i q 4- t~ T ^b T 

for q -- (1 — 2//3)" 1 and r = /?/ 2, we have 

(si 1 p» + ) , - 2/ * 3 ||(u+) 2/ 0|| n 6 < C3upv + 4. £ i“^/*||(„+)2/h||/J/2 V£ > 0 

’ 7l\0 

By taking t = ——j and plugging in our inequality on v, We obtain : 

2C\ a 


28 


sup v < (l/2)supu + + (i/2) 1 “ /3/2 (c 1 ^ 2 |j(u + ) 2 ^|| ;5/ ? (79) 

B n ' 


+( cM )\ i/iu- 

We want to prove the theorem for all p > 0. We will treat separately 
the cases p < n and p> n. 

If p < n, we set fi = 2nfp. In this case we have 

\\(u + ) 2/0 \ \% - l!(« + )IU. 

Plugging this in our inequality on v , we obtain : 


(1/2) sup v < (l/2) , - /J/2 (c 1 dy /2 ||(u + )|| ff + (cid/A)||/|| 6 . 

B 

Consequently, we obtain the following inequality ; 

SUp V < C2 «/*< 

On the sphere i?i/2(0), the cut off function satisfies 


It follows, then 


sup u < sup {vjrj) 

Bj/ 2 (0) Bip(0) 

< 2' ? supv. 

6 

Finally we end up at the desired estimate 

sup it < c 3 {( / (ft 4 ) p ) 1/p + (J/2X)II/H $). 

B,p( o) 

for tt in H ,2,, ‘(fi)nC 2 (fl). The constant C3 above depend only on n,0, -7 and 
6, but is independent of c. 
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On the other ha.nd if p > n, we have : 


2 n /0 < p, V/? > 2. 

Then, it follows (by assuming 0 > 2) 

|B|->/(=»/«|| (ll+) || 2n/oS £ lfil-1/.u.+n^. 

But 


ll« + iij.//» = li(« + ) 2/i, lC5‘ 

and therefore, by processing as before, we obtain the desired estimate 
sup tt < c 4 {( [ (u+y) 1 /” 4- (d/2A)||/|| n jj} 

for m in W 2,n (ft) n <7 2 (ft). The constant c 4 above depends only on n,/3,7 
and S y but is independent of c. 

Transformation & — ♦ x> 

By construction, D tJ = R~ 2 D ij , thus A == R ~ 2 A and 8 = A/£ 2 . In 
addition, we have j.8| = w„(2.R) n and = /T’' i/p |$! PiS . 

Written in term of x, the last inequality becomes 

Z, “ £ c<{( w + (^)ii/iu). 

with c 4 a function of n, 7, <5 ~ <5fi! 2 and p. This is the desired estimate for 
it € W 2 -"(fi) n 
Step 2: 

Now, let u € By density, let, (tt m ) be a sequence of functions 

of C 2 ( 8 ), converging towards u in W 2 ' n (B). The injection of W 2 - n (B) in 
C e ’{B) is continuous, consequently (14,, ) converges uniformly towards u in 
B. We have 


BVffn “ ~ ti J 4 Xit 

> / 4 b(u m - tt). 
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By getting, f m - L(u m - u), we observe by construction that f m con- 
verges towards 0 in L n (J2). As u m £ W 2,n (ft)0 C 2 (fi) and f m — f 4- / m is 
in £ n (fl), the estimate (3) is valid also for u m , so that we have 

sup Um < ^e{(~ f (u^Y) 1/p + “||/|U.b}- (80) 

B R (y) l^i JB * 

Using previous; results and taking the limit, we have : 

sup u<cte{(JL f (u+Y) lh + yii/IU.s}. 

Bp{ y) 1-0 1 JB A 

■ 

Observe also that by replacing u by —u, the theorem, can be extended 
easily to the case c f supersolutions and solutions of the equation : 

Lu = /. 
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